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Abstract. Ultrasound speckle is a granular texture pattern appearing in ultrasound imaging. 
It can be used to distinguish tissues and identify pathologies. Lorentz force electrical 
impedance tomography is an ultrasound-based medical imaging technique of the tissue 
electrical conductivity. It is based on the application of an ultrasound wave in a medium placed 
in a magnetic field and on the measurement of the induced electric current due to Lorentz 
force. Similarly to ultrasound imaging, we hypothesized that a speckle could be observed 
with Lorentz force electrical impedance tomography imaging. In this study, we first assessed 
the theoretical similarity between the measured signals in Lorentz force electrical impedance 
tomography and in ultrasound imaging modalities. We then compared experimentally the 
signal measured in both methods using an acoustic and electrical impedance interface. Finally, 
a bovine muscle sample was imaged using the two methods. Similar speckle patterns were 
observed. This indicates the existence of an “acousto-electrical speckle" in the Lorentz force 
electrical impedance tomography with spatial characteristics driven by the acoustic parameters 
but due to electrical impedance inhomogeneities instead of acoustic ones as is the case of 
ultrasound imaging. 
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Ultrasound (US) B-mode is a widespread medical imaging technique of the tissue acoustical 
impedance. With this technique, an ultrasound wave is transmitted by an acoustic transducer 
and is backscattered by the change in acoustic impedance at tissue interfaces. Measuring the 
analytic signal obtained from ultrasound waves received on the transducer surface leads to the 
reconstruction of these acoustic impedance interfaces. Moreover, when ultrasound waves are 
scattered by acoustic impedance inhomogeneities, their contributions add up constructively 
or destructively on the transducer surface. The total intensity thus varies randomly due to 
these interferences HCobbold, 2007] . The phenomenon appears on reconstructed images as a 
granular pattern with variable size and intensity, which looks spatially random: it is called 
“acoustic speckle” [ [Abbott and Thurstone, 1979] , 

The speckle can be treated as noise and consequently be reduced, or on the 
contrary be considered as a feature [ [Noble and Boukerroui, 20061 ■ Under the first point 
of view, quantitative ultrasound methods were proposed based on the backscattering 
coefficient JLizzi et al., 1983[ |Lizzi et al., 1987[ |Insana et al., 19901 , where the speckle 
is removed to obtain a measure that depends only on tissue acoustic properties; 
see UGoshal et al., 20131 for further references. Under the second point of view, 
quantitative ultrasound methods were developed by exploiting the statistics of the 
echo envelope [Burckhard t, 1978[ |Wagner et al., 19 83] that depend on the speckle and 
tissue characteristics to discriminate tissue types; see [Destrempes and Cloutier, 2010 


Destrempes and Cloutic r72013j Yamaguchi, 2013| for further references. In either case, 


quantitative ultrasound can be used to diagnose pathologies with ultrasound images. 
Another field of application is speckle-tracking that aims at estimating local displacements 
and tissue deformations; see for instance [Ophir et al., 1991, Hein and O’Bri en, 1993[ 


Lubinski et al., 1999f 


On the other hand, the Lorentz force electrical impedance tomography (LFEIT) 
method [ |Montalibet et al., 2001[ |Grasland-Mongrain et al., 2013 1, also known as magneto- 
acousto electrical tomography [[Haider et al., 2008ft , is a medical imaging technique producing 
electrical conductivity images of tissues flWen et al., 1998[|Roth, 2Qlll|Ammari et al., 2014[ [. 
With this technique, an ultrasound wave is transmitted by an acoustic transducer in a 
biological tissue placed in a magnetic field. The movement of the tissue in the magnetic 
field due to the ultrasound propagation induces an electric current due to Lorentz force. The 
measurement of this current with electrodes in contact with the sample allows to reconstruct 
images of electrical impedance interfaces. 

The goal of this work was to assess the presence of an “acousto-electrical” speckle 
in the LFEIT technique, similar to the acoustic speckle in US imaging, as suggested in 
a previous work [ [Haider et al., 2008| |. The first part of this work presents the theoretical 
similarity of measured signals in these two imaging techniques. This similarity is then 
observed experimentally on an acoustic and electrical interface. Finally, a bovine sample 
is imaged using both methods to observe the two types of speckle. 
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2. Theoretical background 

The goal of this section is to compare the mathematical framework of the LFEIT method with 
that describing radio-frequency signals in US imaging. 

2.1. Local current density in Lorentz force electrical impedance tomography 



Magnetic 

Field 


Figure 1 . An ultrasound wave is transmitted in a conductive medium placed in a magnetic 
field. This induces an electric current I to t due to Lorentz force. This current separates 
in two components, Ii ns which stays inside the medium, and I meas which is measured 
by electrodes in contact. The measured signal allows to reconstruct images of electrical 
impedance interfaces. 


The LFEIT principle is illustrated in Figure |T] In this technique, an ultrasound wave is 
transmitted along a unit vector e z in a conductive medium placed in a magnetic field along 
e x of intensity B x . This induces an electric current due to Lorentz force in the direction 
e y = e z x e x . We assume that the conductive medium is moved by the ultrasound wave along 
the e z direction so that all particles within the medium are moved with a mean velocity of 
amplitude v z parallel to e 2 [|Mari et al., 2009|| . Since the medium is placed in a magnetic field 
B x , a particle k of charge q k is deviated by a Lorentz force F k = q k v z B x e y . Using Newton’s 
second law, the velocity u k of the charged particle can be calculated as: 

Ufc t] Z G Z T p, k v z B x & y ^ (1) 

where p, k is the mobility of the particle /;:. 

The density of current j, defined as q k per unit volume, is then equal to the sum 
of ffk QkVz^z and ff k q k p k v z B x e y per unit volume. If the medium is assumed electrically 
neutral, i.e. q k = 0, the first term is equal to zero. By introducing the electrical 
conductivity o, defined as q k p k per unit volume, the second term is equal to ov z B x e y . 
Equation (JTJ) can consequently be written as: 


j = crv z B x e y . 


( 2 ) 
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The local density of current j is thus proportional to the electrical conductivity of the sample 
a, to the component B x of the applied magnetic field, and to the local sample speed v 


[Montalibet et al., 2001]. 


2.2. Measured signal in Lorentz force electrical impedance tomography 

Equation ([2]) is however local. We consider, as a first approximation, the ultrasound beam as a 
plane wave along e z inside a disk of diameter W, with the velocity v z = 0 outside the beam. 
The total current I tot induced by the Lorentz force is then equal to the average of the current 
density j computed over all surfaces S inside the ultrasound beam that are perpendicular to 


[ |Montalibet, 2002| : 


~ w 



j.dS dy = w 



av z B x dz dx dy. 


(3) 


As electrodes are placed outside the beam, we assume they measure only a fraction a of the 
total current I tot induced by the Lorentz force. The current I tot then induces a current outside 
the beam which separates in two components: a current I ms that stays inside the medium 
and operates through a resistance R, ms , and a current I mea s that is measured by electrodes 
and operates through a resistance R meas - Due to Ohm’s law, which states that L m Jl nis = 
ImeasRmeasi the Coefficient (\ , defined as Irneas/ftoti is also equal to Rins/ (Rmeas Rins) • 
This coefficient, rather difficult to estimate, depends consequently on a few parameters: the 
size and location of electrodes, and the circuit and medium electrical impedance. We can 
nevertheless conclude that the lower the resistance R meas compared to Ri ns , the higher a is, 
and the higher is the measured electrical current. 

Lor a linear and inviscid medium, the medium velocity v z due to ultrasound wave 
propagation is related to the ultrasound pressure p and the medium density p (assumed here 
constant over time but not necessarily over space) using the identity dVz g t ’ z) = ~y 9p j)f' > 
(having assumed a plane wave). Lrom these considerations on the coefficient a and the 
medium velocity v z , the measured current I meas is equal to: 

rZ2 , rt 


, , a 

•« = w 



crBy 


1 dp(r, z ) 
dz 


dr ) dz dx dy , 


(4) 


/ X Oi 

= w 



' 21 'J — OO P 

where z\ and z^ are the boundaries of the studied medium along e z . 

Considering a progressive acoustic wave that propagates only along the direction e z and 
is not attenuated in the measurement volume, we may assume that p(r, z) is of the form 
P(t — z/c) where c is the speed of sound in the medium. Equation [4] can consequently be 
written as: 

1 

- 2 ! PC 

where P(— oo) = 0 since the transmitted pulse is of finite length. 

Assuming B x and c constant over space (c varies typically from -5 to +5% between soft 
biological tissues [|Hill et al., 2004|| ), and replacing the integration over z by an integration 
over t = z/c, equation [5] becomes: 

B roo 

Imeas (f) = / H z (j)P{t - t) dr , ( 6 ) 


z/c) — P(—oo)^j dz dx dy , 


(5) 
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where H z (r) = f f H(x', y', z' = CT)dx'dy', H is equal to - p within the studied medium and 
0 elsewhere. Since the medium density vary typically by a few percent between different 


soft tissues [[Cobbold, 20071 , while the electrical conductivity of soft tissues can vary up to 
a few tens [ [Gabriel et al., 19961 , H can be seen mostly as a variable describing the electrical 
conductivity of the tissue. One can show that if the DC component of the transmitted pressure 
is null, then I meas is null whenever the electrical conductivity is constant [ [Wen et al., 1998Q . 
Thus, equation [6] shows that the measured electric current is proportional to the convolution 
product of a function H z of the electrical conductivity with the axial transmitted pressure 
wave P. 


2.3. Measured signal in B-mode ultrasound imaging 

Ultrasound imaging is based on the measurement of reflections of the transmitted acoustic 
wave (i.e., backscattering, diffraction or specular reflection). We assume that acoustic 
inhomogeneities are scattering only a small part of the transmitted pressure, so that scattered 
waves have a negligible amplitude compared to the main acoustic beam, and that the 
diffraction and attenuation are small in the region of interest. The voltage RF(t ) measured at 
time t on the transducer, known as the radiofrequency signal, is then equal to: 


/ OO 

T z (r)P(t 

-OO 


r) dr, 


(7) 


where D is a constant of the transducer related to the acousto-electric transfer function, 
T,(t) = f J T{x ', y', z' = ct/2) dx' dy' with T the continuous spatial distribution of point 
scatterers, P(t) is the axial transmitted pressure wave (the pulse shape), and r = 2 zjc 
(the factor 2 takes into account round-trip wave propagation in ultrasound imaging). In the 
latter expression of T z (r), the double integral is performed over a disk of diameter W (the 
ultrasound beam is considered as a plane wave along e z inside a limited width W as in the 


previous section) [ jBamber and Dickinson, 19801 - The envelope of RF(t ) represents the A- 
mode ultrasound signal. Note that under a more realistic incident pressure wave model, one 
can consider pulse beam profiles in the definition of H z (at emission) and T z (at emission 
and at reception), which can be deduced by multiplying the pulse shape with the incident 
pulse beam profile. Also, attenuation can be taken into account by convolution of Eq. [7] (and 
similarly, of Eq. [6]) with an attenuation function. 

Hence, the received signal in US imaging is proportional to the convolution product of 
the axial spatial distribution of point scatterers T z with the transmitted pulse shape P. If 
scatterers are small compared to the ultrasound wavelength, reflected waves can interfere. 
When the envelope of the signal is calculated and a B-mode image is formed line by 
line, this phenomenon appears on the image as a granular texture called acoustic speckle 


[Abbott and Thurstone, 19791. 


Apart from a proportionality coefficient, equations ([6]) and ([7]) have a similar form, where 
T and H play an analogous role. A speckle phenomenon is thus expected to be observed in 
the LFEIT technique. Its nature would be electrical, because mainly related to electrical 











Acousto-electrical speckle pattern in Lorentz force electrical impedance tomography 6 

conductivity inhomogeneities, and also acoustic, because spatial characteristics are related to 
the acoustic wavelength. 

3. Materials and methods 

Two experiments were performed in this study. The first experiment aimed at comparing a 
LFEIT signal and a US signal on a simple acoustic and electrical conductivity interface using 
the same acoustic transducer. This approach was used to perform a test with a large change 
in both electrical conductivity and acoustic impedance. Thus, tissue functions T z and H z 
could be considered as square pulse functions with strong gradients at the interface locations. 
According to the above hypotheses, the two scan line signals are expected to be identical apart 
from a proportionality coefficient in this special case. 

The goal of the second experiment was to observe a complex biological tissue with the 
two imaging techniques and using the same acoustic transducer. A speckle pattern of similar 
spatial characteristics is expected to be observed on both types of images, since they are 
related to the acoustic wavelength and beam width, but different, because the nature of the 
imaged parameter is acoustic in one case and electrical in the other case. 

3.1. Measured signal on an acousto-electrical interface 

The experiment setup is illustrated in Figure [2j A generator (HP33120A, Agilent, Santa 
Clara, CA, USA) was used to create 0.5 MHz, 3 cycles sinusoid bursts at a pulse repetition 
frequency of 100 Hz. This excitation was amplified by a 200 W linear power amplifier (200W 
LA200H, Kalmus Engineering, Rock Hill, SC, USA) and sent to a 0.5 MHz, 50 mm in 
diameter transducer focused at 210 mm and placed in a degassed water tank. The peak- 
to-peak pressure at the focal point was equal to 3 MPa. A 4x15x20 cm 3 mineral oil tank 
was located from 15 to 35 cm away from the transducer in the ultrasound beam axis. This 
oil tank was used to decrease the loss of current from the sample to the surrounding medium 
and to increase consequently the signal-to-noise ratio, but was not mandatory. It was inserted 
in a 300±50 mT magnetic field created by a U-shaped permanent magnet, composed of two 
poles made of two 3x5x5 cm 3 NdFeB magnets (BLS Magnet, Villers la Montagne, France) 
separated by a distance of 4.5 cm. 

The tested medium was a 5x10x10 cm 3 10% gelatin filled with 5% salt sample placed 
inside the oil tank from 30 to 40 cm away from the transducer, and presented consequently a 
strong acoustic and electrical interface. A pair of 3x0.1x10 cm 3 copper electrodes was placed 
in contact with the sample, above and under it, respectively. The electrodes were linked 
through an electrical wire to a 1 MV/A current amplifier (HCA-2M-1M, Laser Components, 
Olching, Germany). A voltage amplifier could also be used, but the current amplifier presents 
a smaller input impedance (a few Ohms vs 50 Ohms) and increases consequently the amount 
of current measured by the electrode, as depicted by the factor a in equation |4j The 
signals were then measured by an oscilloscope with 50 U input impedance (WaveSurfer 422, 
LeCroy, Chestnut Ridge, NY, USA) and averaged over 1000 acquisitions. US signals were 
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Figure 2. A transducer is transmitting ultrasound pulses toward a sample placed in an oil tank 
placed in a magnetic field. The induced electric current is received by two electrodes in contact 
respectively with two sides of the gelatin. 


To quantify similarities between the two signals, we computed the correlation coefficient 
between them. Such coefficient is equal to 1 when the two signals are directly proportional 
and 0 when they are uncorrelated. 


3.2. Observation of speckle in a biological tissue 

The same apparatus as in the first experiment was used, but the gelatin sample was replaced 
by a 2x6x6 cm 3 piece of bovine rib muscle purchased at a grocery store. It presented many fat 
inclusions, as pictured in Figure [3] B-mode images were produced line by line by moving the 
transducer along the e y direction by 96 steps of 0.5 mm. Acoustic and electrical signals were 
post-processed using the Matlab software (The Math Works, Natick, MA, USA) by computing 
the magnitude of the Hilbert transform of the signal [|Roden, 19911 , and were displayed with 
grayscale and jet color images, respectively. 


4. Results 


Figure |4]-(A) presents the electrical signal measured by electrodes from the first phantom 
interface, 195 to 220 microseconds after acoustic transmission, corresponding to a distance 
of 30 cm. Figure |4]-(B) depicts the signal acquired by the US transducer from the phantom 
interface, 395 to 420 microseconds after acoustic transmission, corresponding to a distance 
of 30 cm (back and forth). Both signals consist in three to four cycles at a central frequency 
of 500 kHz. The correlation coefficient between both signals is 0.9297, indicating a high 
similarity. This similarity was observed at the second phantom interface, with a correlation 
coefficient between both signals equal to 0.9179. LFEIT and US images of the phantom can 


be seen in [Grasland-Mongrain et al., 2013]. 
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Figure 3. Picture of the 2x6x6 cm 3 imaged bovine rib muscle sample, with many fat 
inclusions. 



Figure 4. (A) Electrical signal acquired by electrodes from the first phantom interface, 195 to 
220 ps after acoustic wave transmission. The signal is made of three to four cycles at 500 kHz. 
(B) Electrical signal acquired by the acoustic transducer from the first phantom interface, 395 
to 420 ps after acoustic wave transmission. The signal is also made of three to four cycles at 
500 kHz. 


Figures [5]- (A) and -(B) present the Lorentz force electrical impedance tomography image 
and the ultrasound image, respectively, of the bovine muscle sample. The amplitude varied 
from -2 to -2.8 dB with the first technique and from -1.5 to -2.5 dB with the other (0 dB being 
a measured amplitude of 1 V). Main interfaces of the medium can be retrieved on both images, 
as previously shown [Grasland-Mongrain et al ., 2013[ , even if signals at the boundaries and 
the interior of the bovine sample were of similar amplitude. A speckle pattern was present 
in both images. The typical spots were of size of same order of magnitude; i.e., 5-8 mm in 
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the Z-direction and 10-18 mm in the Y-direction, but their spatial distributions were different, 
as expected due to the difference between the electrical and acoustic inhomogeneities. For 
each image, the signal-to-noise ratio was estimated as the base 10 logarithm of the ratio of 
the mean square amplitude of the RF signals in the 22.5-26 cm zone with the mean square 
amplitude in the 19-20 cm noisy zone. The signal-to-noise ratio was lower by 0.9 dB in the 
LFEIT image compared to the US image. 



z (cm) z (cm) 

Figure 5. (A) Lorentz force electrical impedance tomography (LFEIT) image of the bovine 
muscle sample. (B) Ultrasound (US) image of the same bovine rib muscle sample. A speckle 
pattern can be seen inside the medium. 


5. Discussion 


The gelatin phantom used in the first experiment presented an interface of acoustic and 
electrical impedances. According to the correlation coefficient, the measured signals were 
very similar, which is a first indication of the validity of the approach presented in the 
theoretical section: the reflected wave is proportional to the convolution product of the 
acoustic impedance distribution with the transmitted ultrasound pulse shape, while the 
induced electric current is proportional to the convolution product of the electrical impedance 
distribution with the transmitted ultrasound pulse shape. 

The second experiment showed two images of granular pattern. This pattern does not 
represent macroscopic variations of acoustic or electrical impedances, and we interpreted it 
as speckle. The granular pattern appeared visually with similar characteristics of size and 
shape in both images. Note that the size of the speckle spots along the ultrasound beam was 
different than in the orthogonal direction, because the first is mainly related to the acoustic 


wavelength and the second to the ultrasound beam width [[O’ Brien , 1993 j. The spots were 
quite larger than those usually seen in clinical ultrasound images, due to the characteristics of 
the instrument (three cycles at 500 kHz, 1.5 cm beam width). 

Although the bright spots were of similar size, their precise locations in the two images 
were different. This shows that the observed speckle reveals information of different nature 
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in the two modalities: acoustic or electrical inhomogeneities. In this case, the bovine meat 
sample presented not only large layers of fat but also small inclusions of adipose tissues 
whose electrical conductivity differs from muscle (differences can be ten time higher at 500 
kHz [ [Gabriel et al., 19961 ), whereas they have a close acoustic impedance (differences smaller 
than 10% [jCobbold, 20071 ). 

We introduce henceforth the term “acousto-electrical speckle”. It is justified by the fact 
that its spatial characteristics are related to the acoustic parameters, especially the ultrasound 
wavelength, while its nature is related to the electrical impedance variation distribution H. 
The existence of this speckle could allow using speckle-based ultrasound techniques in the 
Lorentz force electrical impedance tomography technique, for example compound imaging, 
speckle-tracking algorithm or quantitative ultrasound for tissue characterization purposes 
HO’ Donnell et al., 1994[ |Jespersen et al., 1998[ [Mamou and Oelze, 2013] , These techniques 
have however not be applied in this study because of the low spatial resolution of the images 
due to the low frequency transducer used. 

This study shows that the electrical impedance inhomogeneities can be studied using 
LFEIT at a scale controlled by the acoustic wavelength instead of the electromagnetic 
wavelength, which is 5 orders of magnitude larger at a same frequency, and would be 
prohibitive in the context of biological tissues imaging. 

These inhomogeneities should also be observed in a “reverse” mode (terminology 
from Wen et al. []Wen et al., 1998| ) called Magneto-Acoustic Tomography with Magnetic 
Induction, where an electrical current and a magnetic field are combined to produce ultrasound 
waves [Rot h et al., 1994[ |Xu and He, 2005[ |Ammari et al., 20091 . However, presence of 
speckle in this last technique has not yet been demonstrated and which of the two methods 
will be most useful in biomedical imaging is not clear [ [Roth, 2011 1. 


6. Conclusion 

In this study, the similarity between two imaging modalities, the Lorentz force electrical 
impedance tomography and ultrasound imaging, was assessed theoretically. This similarity 
was then observed experimentally on a basic acoustic and electrical interface with both 
methods. Then, the two techniques were used to image a biological tissue presenting 
many acoustic and electrical impedance inhomogeneities. The speckle pattern formed in 
both images exhibited similar spatial characteristics. This suggests the existence of an 
“acousto-electrical speckle” with spatial characteristics driven by acoustic parameters but 
due to electrical impedance variation distribution. This allows considering the use of 
ultrasound speckle-based image processing techniques on Lorentz force electrical impedance 
tomography data and to study electrical inhomogeneity structures at ultrasound wavelength 
scale. 
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